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ABSTRACT 


Many algorithms used for data smoothing, data classification and 
error detection require the calculation of the distance from a 
point to the polynomial interpolating its 2n neighbors (n on 
each side) as is clearly the case in [1]. This computation, if 
performed naively, would require the solution of a system of equations 
and could create numerical problems. It is the purpose of this note 
to show that if the data is equally spaced, then this calculation 
can be performed using a simple recursion formula. 
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A GEOMETRICAL INTERPRETATION OF THE 2n-th CENTRAL DIFFERENCE* 


R. A. Tapia** 

1. Introduction . Many algorithms used for data smoothing, data classifi- 
cation and error detection require the calculation of the distance from a 
point to the polynomial interpolating its 2n neighbors (n on each side) as 
is clearly the case in [1]. This computation, if performed naively, would 
require the solution of a system of equations and could create numerical 
problems. It is the purpose of this note to show that if the data is equally 
spaced, then this calculation can be performed using a simple recursion 
formula. 


2. The 2n-th Central Difference . Given (wj : i = 0, +1, +2 , . . . jcR 1 and a 

scalar h we may define the second difference operator 

w 1+1 - 2w t + w t _i 
D (wj ;h) = jp 


Consider Y = ( (x t ,yj )c: R 2 : i = 0,+l,+2, — } with x^ +1 - x t = h for all i 
Recall that by the 2n-th central difference of Y at the i-th point we 
mean 


r ( i 2n) = D(rP n ~ 2) ;h), n = 1,2,3,... 


with r$ 0 ^ = y l for all i . 

By the 2n-th canonical central difference r^ 2n ^ we mean the 2n-th 
central difference of the set 


C Cjh, 6 0 J ) : j = 0,+l,+2,...J 
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at (0,1) (6j j is the Kronecker delta). 

A well-known property of the 2n-th central difference is that it is 
a good approximation to the 2n-th derivative for small h, i.e., if 
Yj = f (x s ) and f is 2n times differentiable, then rj 2n ^ -> f^ 2n ^(xj ) as 
h 0 . We offer the following alternate and very useful interpretation of 
the 2n-th central difference. 


Theorem . Let P(x) be the unique polynomial of degree 2n - 1 interpolating 
the points 

((*i+j J = ± 1 >±2, • • • ,±n] (n > 1) . 


Then 


r l 

Yi “ ? ( x i ) = r ' ( a n ) for any h $ 0 . 
Proof . We merely sketch a proof. Consider the following table. 


n 

Yi-4 

yi-s 

yi-2 

TABLE 1 
yi-1 Vi 

yi+i 

yi+2 

yi+a 

yt+4 

0 

0 

0 

0 

0 1 

0 

0 

0 

0 

1 

0 

0 

0 

1 -2 

1 

0 

0 

0 

2 

0 

0 

1 

-4 6 

-4 

1 

0 

0 

3 

0 

1 

-6 

15 -20 

15 

-6 

i 

0 

4 

1 

-8 

28 

-56 70 

-56 

28 

-8 

1 

Table 

1 is constructed 

by taking for the 

n-th 

row the 

second 

difference > 

the 

(n-l)-th 

row with 

h = 

1, for n = 

1,2,. 

• • • Let a n j 

denote the 


element in Table 1 in the n-th row and in the column corresponding to yj . 
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From Table 1 we obtain r[ zn) = £ a nj y 3 • Moreover r (2r ° = a^ t . 

1 

We now construct Table 2 from Table 1 by dividing by the column 
corresponding to y i , i.e., by dividing the n-th row by ^ r^ 2 * 1 -* . 

TABLE 2 


n 

y t-4 Yi-3 

yi-a 

Yi-l 

yi 

yi+i 

yi + 2 

y i+ 3 

y t +4 

0 

0 

0 

0 

0 

i 

0 

0 

0 

0 

1 

0 

0 

o 

-1/2 

i 

-1/2 

0 

0 

0 

2 

0 

0 

1/6 

-2/3 

i 

-2/3 

1/6 

0 

0 

3 

0 

-1/20 

3/10 

-3/4 

i 

-3/4 

3/10 

-1/20 

0 

4 

1/70 

-4/35 

4/10 

-4/5 

i 

-4/5 

4/10 

-4/35 

1/70 

As in 

Table 1 let 

b n , 

denote 

a typical 

el ement 

of Table 2. 

It can be shown 

that 

b n,l+j for 

-n < 

j < n 

is the 

value at 

x i of 

the Lagrangian poly- 


nomial of degree 2n - 1 which interpolates the n points on each side of 
x r and is equal to 1 at x l+3 and 0 at the remaining points. Once 
this has been established, using Lagrange's formula it follows that 

y l - P(xj) -Hjyj ; 

E b nj . This proves the theorem . 
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